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Abstract 

The spectra of strange hadrons are measured in proton-proton collisions, recorded by 
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The Kg, A, and particles and their antiparticles are reconstructed from their decay 
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by up to a factor of three. 
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1 Introduction 

Measurements of particle yields and spectra are an essential step in understanding proton- 
proton collisions at the Large Hadron Collider (LHC). The Compact Muon Solenoid (CMS) 
Collaboration has published results on spectra of charged particles at centre-of-mass ener- 
gies of 0.9, 2.36, and 7 TeV[l, 2J. In this analysis the measurement is extended to strange 
mesons and baryons (Kg, A, E^)[^at centre-of-mass energies of 0.9 and 7 TeV. The investi- 
gation of strange hadron production is an important ingredient in understanding the nature 
of the strong force. The LHC experiments ALICE and LHCb have recently reported results 
on strange hadron production at a/s = 0.9 TeV ||3l|4l. In addition to results at ^/s = 0.9 TeV, 
we also present results at = 7 TeV, opening up a new energy regime in which to study 
the strong interaction. As the strange quark is heavier than up and down quarks, production 
of strange hadrons is generally suppressed relative to hadrons containing only up and down 
quarks. The amount of strangeness suppression is an important component in Monte Carlo 
(MC) models such as pythia |5J and hijing/bb |6]. Because the threshold for strange quark 
production in a quark-gluon plasma is much smaller than in a hadron gas, an enhancement 
in strange particle production has frequently been suggested as an indication of quark-gluon 
plasma formation [7J. This effect would be further enhanced in baryons with multiple strange 
quarks. While a quark-gluon plasma is more likely to be found in collisions of heavy nuclei, 
the enhancement of strange quark production in high energy pp collisions would be a sign of a 
collective effect, according to some models ||81|9J. In contrast, recent Regge-theory calculations 
indicate little change in the ratio of Kg to charge particle production with increasing collision 
energy IITOllTTl . Thus, these measurements can be used to constrain theories, provide input for 
tuning of Monte Carlo models, and serve as a reference for the interpretation of strangeness 
production results in heavy-ion collisions. 

Minimum bias collisions at the LHC can be classified as elastic scattering, inelastic single- 
diffractive dissociation (SD), inelastic double-diffractive dissociation, and inelastic non-diffractive 
scattering. The results presented here are normalized to the sum of double-diffractive and 
non-diffractive interactions, referred to as non-single-diffractive (NSD) interactions III] |2]. This 
choice is made to most closely match the event selection and to compare with previous experi- 
ments, which often used similar criteria. The Kg, A, and are long-lived particles (ct > 1 cm) 
and can be identified from their decay products originating from a displaced vertex. The par- 
ticles are reconstructed from their decays: Kg — > tt+tt^, A — > pn^, and — > An^ over the 
rapidity range \y\ < 2, where the rapidity is defined as y = ^ In f^^/ £ is the particle energy, 
and Pl is the particle momentum along the anticlockwise beam direction. For each particle 
species, we measure the production rate versus rapidity and transverse momentum pj, the 
average pj, the central production rate ^\yrvO> and the integrated yield for \y\ < 2 per NSD 
event. We compare our measurements to results from Monte Carlo models and lower energy 
data. 

2 CMS experiment and collected data 

CMS is a general purpose experiment at the LHC iT2\ . The silicon tracker, lead-tungstate crystal 
electromagnetic calorimeter, and brass-scintillator hadron calorimeter are all immersed in a 
3.8 T axial magnetic field while muon detectors are interspersed with flux return steel outside of 
the 6 m diameter superconducting solenoid. The silicon tracker is used to reconstruct charged 
particle trajectories with < 2.5, where the pseudorapidity is defined as rj = — In tan 6 

^Particle-conjugate states are implied throughout this paper. 
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3 Strange particle reconstruction 



being the polar angle with respect to the anticlockwise beam. The tracker consists of layers 
of 100 X 150 }im^ pixel sensors at radii less than 15 cm and layers of strip sensors, with pitch 
ranging from 80 to 183 }im, covering radii from 25 to 110 cm. In addition to barrel and endcap 
detectors, CMS has extensive forward calorimetry including a steel and quartz-fibre hadron 
calorimeter (HF), which covers 2.9 < \rj\ < 5.2. The data presented in this paper were collected 
by the CMS experiment in spring 2010 from proton-proton collisions at centre-of-mass energies 
of 0.9 and 7TeV during a period in which the probability for two collisions in the same bunch 
crossing was negligible and the bunch crossings were well separated. 

The online selection of events required activity in the beam scintillator counters at 3.23 < < 
4.65 in coincidence with colliding proton bunches. The offline selection required deposits of 
at least 3 GeV of energy in each end of the HF preferentially selecting NSD events. A pri- 
mary vertex reconstructed in the tracker was required and beam-halo and other beam-related 
background events were rejected as described in Ref. [IJ. The data selected with these cri- 
teria contain 9.08 and 23.86 million events at 0.9 and 7 TeV, corresponding to approximate 
integrated luminosities of 240 and 480 f/b^^, respectively. To determine the acceptance and ef- 
ficiency, minimum-bias Monte Carlo samples were generated at both centre-of-mass energies 
using PYTHIA 6.422 |5| with tune D6T [13 J. These events were passed through a CMS detector 
simulation package based on GEANT 4 ||T4[. 

3 Strange particle reconstruction 

Ionization deposits recorded by the silicon tracker are used to reconstruct tracks. To maxi- 
mize reconstruction efficiency, we use a combined track collection formed from merging tracks 
found with the standard tracking described in Ref. [l5j and the minimum bias tracking de- 
scribed in Ref. [IJ. Both tracking collections use the same basic algorithm; the differences are in 
the requirements for seeding, propagating, and filtering tracks. 

As described in Ref. [15J, the Kg and A (generically referred to as V*^) reconstruction combines 
pairs of oppositely charged tracks; if the normalized X'^ of the fit to a common vertex is less than 
7, the candidate is kept. The primary vertex is refit for each candidate, removing the two tracks 
associated with the candidate. The next two paragraphs describe the selection of candidates 
for measurement of V*^ and properties, respectively. Selection variables are measured in 
units of c, the calculated uncertainty including all correlations. 

To remove Kg particles misidentified as A particles and vice versa, the Kg (A) candidates must 
have a corresponding pn^ {n^n^) mass more than 2.5a away from the world-average A(Kg) 
mass. The production cross sections we measure are intended to represent the prompt produc- 
tion of Kg and A, including strong and electromagnetic decays. However, particles can also 
be produced from weak decays and from secondary nuclear interactions. These unwanted con- 
tributions are reduced by requiring that the V*^ momentum vector points back to the primary 
vertex. This is done by requiring the 3D distance of closest approach of the V*^ to the primary 
vertex to be less than 3a. To remove generic prompt backgrounds, the 3D V*^ vertex separation 
from the primary vertex must be greater than 5a and both V*^ daughter tracks must have a 3D 
distance of closest approach to the primary vertex greater than 3a. With the above selection, 
the background level for low transverse-momentum A candidates remains high. Therefore, 
additional cuts are applied to A candidates with pj < 0.6 GeV/ c: 

• 3D separation between the primary and A vertices > 10a (instead of > 5a), 

• transverse (2D) separation between the pp collision region (beamspot) and A vertex 
> 10a (instead of no cut), where the uncertainty is dominated by the A vertex, and 
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• 3D impact parameter of the pion and proton tracks with respect to the primary ver- 
tex > (7 — 2|i/|)cr (instead of > 3cr) where y is the rapidity of the A candidate. The 
rapidity dependence is a consequence of the observation that, for the low transverse 
momentum candidates, large backgrounds dominate at small rapidity, while low 
efficiency characterizes the large rapidity behaviour. 

The resulting mass distributions of Kg and A candidates from the 0.9 and 7 TeV data are shown 
m Figs. |l| and g The 7T+ n mass distribution is fit with a double Gaussian (with a common 
mean) signal function plus a quadratic background. The pn^ mass distribution is fit with 
a double Gaussian (common mean) signal function and a background function of the form 
Aq^, where q = Mpjj^- — {nip + nij^--), Mpj^- is the pn^ invariant mass, and A and B are free 
parameters. The fitted Kg (A) yields at ^/s = 0.9 and 7 TeV are 1.4 x 10^ (2.8 x 10^) and 6.5 x 
lO'^ (1.5 X 10^), respectively 




Figure 1: The n^rc^ invariant mass distributions from data collected at = 0.9 TeV (left) and 
7 TeV (right). The solid curves are fits to a double Gaussian and quadratic polynomial. The 
dashed curves show the quadratic background contribution. 



To reconstruct the S^, charged tracks of the correct sign are combined with A candidates. The 
probability of the fit to a common vertex for the A and the charged track must be greater 
than 5%. In this fit, the A candidate is constrained to have the correct world-average mass IIT6l . 
The primary vertex is refit for each candidate, removing all tracks associated with the H^. 
The candidates must then pass the following selection criteria: 

• 3D impact parameter with respect to the primary vertex > 2a for the proton track 
from the A decay, > 3a for the track from the A decay, and > 4(7" for the tt^ track 
from the decay, 

• invariant mass from the hypothesis for the tracks associated with the A can- 
didate at least 20 MeV/ away from the world-average Kg mass, 

• 3D impact parameter of the candidate with respect to the primary vertex < 3a, 

• 3D separation between A vertex and primary vertex > \0a, and 

• 3D separation between vertex and primary vertex > 2a. 

The mass distributions of candidates from the i/s = 0.9 and 7 TeV data are shown in Fig.js] 
The Att^ mass is fit with a double Gaussian (with a common mean) signal function and a 
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4 Efficiency correction 




pji invariant mass [MeV/c^] p% invariant mass [MeV/c^] 

Figure 2: The pn^ invariant mass distributions from data collected at a/s = 0.9 TeV (left) and 
7TeV (right). The solid curves are fits to a double Gaussian signal and a background func- 
tion given by Aq^, where q = Mpj^- — {nip + nij^-). The dashed curves show the background 
contribution. 



background function of the form Aq^^^ + Bq^^'^, where q = M/^jj^- — (/Wa + mj^--) and Mati- is 
the Att^ invariant mass. The fitted yields at ^^s = 0.9 and 7 TeV are 6.2 x 10^ and 3.4 x 10^, 
respectively. 

4 Efficiency correction 

The efficiency correction is determined from a Monte Carlo simulation which is used to mea- 
sure the effects of acceptance and the efficiency for event selection (including the trigger) and 
particle reconstruction. The Monte Carlo samples are reweighted to match the observed track 
multiplicity in data, as this has been shown to be an important component of the trigger effi- 
ciency |ITJ|2l. This is referred to as track weighting. The efficiency correction also accounts for 
the other decay channels of the strange particles that we do not attempt to reconstruct, such as 

The efficiency is given by the number of reconstructed particles divided by the number of 
generated particles, subject to two modifications. Firstly, the efficiency correction is used to ac- 
count for candidates from SD events. As the results are normalized to NSD events, candidates 
from SD events which pass the event selection must be removed. This is done by defining the 
efficiency as the number of reconstructed candidates in all events divided by the number of 
generated candidates in NSD events. Secondly, the efficiency is modified to account for the 
small contribution of reconstructed non-prompt strange particles which pass the selection cri- 
teria. This is only an issue for the A particles which receive contributions from S and Q decays. 
Since these non-prompt A particles are present in both the MC and data, we modify the effi- 
ciency to remove this contribution by calculating the numerator using all of the reconstructed 
strange particles and the denominator with only the prompt generated strange particles. As 
the MC fails to produce enough E particles (see Section |6]l, the non-prompt A's are weighted 
more than prompt A's in the efficiency calculation. 

The results of this analysis are presented in terms of two kinematic distributions: transverse 
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Figure 3: The Att^ invariant mass distributions from data collected at v^s = 0.9 TeV (left) 7TeV 
(right). The solid curves are fits to a double Gaussian signal and a background function given 
by Aq^^^ + Bq^^'^, where q = M^tt- — (wza + nij^-). The dashed curves show the background 
contribution. 

momentum and rapidity. For all modes, is divided into 10 equal size bins from to 2 and 
px is divided into 20 equal size bins from to 4 GeV/c plus one bin each from 4 to 5 GeV/c and 
5 to 6 GeV/c. In addition, the V° modes also have 6-8 GeV/c and 8-10 GeV/c px bins. All results 
are for particles with \y\ < 2. 

The efficiency correction for the modes uses a two-dimensional binning in pj and \y\. Thus, 
the data are divided into 240 bins in the \y\, px plane. The invariant mass histograms in each 
bin are fit to a double Gaussian signal function (with a common mean) and a background 
function. In bins with few entries, a single Gaussian signal function is used. For the A sample, 
some bins are merged due to sparse populations in \y\,pj space. The merging is performed 
separately when measuring \y\ and px such that the merging occurs across px and \y\ bins, 
respectively. The efficiency from MC is evaluated in each bin and applied to the measured yield 
to obtain the corrected yield. The two-dimensional binning used for the V*^ efficiency correction 
greatly reduces problems arising from remaining differences in production dynamics between 
the data and the simulation. The much smaller sample of candidates prevents the use 
of 2D binning. Thus, the data are divided into \y\ bins to measure the \y\ distribution and 
into Px bins to measure the px distribution. However, the MC spectra do not match the data. 
Therefore, each Monte Carlo particle is weighted in px i\y\) to match the distribution in 
data when measuring the efficiency versus \y\ (px)- Thus, the MC and data distributions are 
forced to match in the variable over which we integrate to determine the efficiency. We refer 
to this as kinematic weighting. The efficiencies for all three particles are shown versus \y\ 
and Px in Fig. [1] The efficiencies (for particles with \y\ < 2) include the acceptance, event 
selection, reconstruction and selection, and also account for other decay channels. The increase 
in efficiency with px is due to the improvement in tracking efficiency as track px increases and 
to the selection criteria designed to remove prompt decays. The slight decrease at high px is 
due to particles decaying too far out to have reconstructed tracks. While there is no centre-of- 
mass energy dependence on the efficiency versus px, particles produced at a/s = 7 TeV have a 
higher average-px, resulting in a higher efficiency when plotted versus rapidity. 

As a check on the ability of the Monte Carlo simulation to reproduce the efficiency, the (well- 
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Figure 4: Total efficiencies, including acceptance, trigger and event selection, reconstruction 
and particle selection, and other decay modes, as a function of \y\ (left) and pj (right) for Kg, 
A, and produced promptly in the range \y\ < 2. Error bars come from MC statistics. 



known) Kg, A, and lifetimes are measured. For the Kg measurement, the data are divided 
into bins of pj and ct, where ct is calculated as = cmL/p where m, L, and p are, respectively, 
the mass, decay length, and momentum of the particle. In each bin the data is corrected by 
the MC efficiency and the corrected yields summed in pj to obtain the ct distribution. Due to 
smaller sample sizes, the A and yields are only measured in bins of ct. Using the kinematic 
weighting technique, the MC efficiency in each bin of ct is calculated with the pj spectrum 
correctly weighted to match data. The corrected lifetime distributions, shown in Fig. |5j display 
exponential behaviour. The vertex separation requirements result in very low efficiencies and 
low yields in the first lifetime bin and are thus expected to have some discrepancies. An actual 
measurement of the lifetime would remove this issue by using the reduced proper time, where 
one measures the lifetime relative to the point at which the particle had a chance to be recon- 
structed. The measured values of the lifetimes are also reasonably consistent with the world 
averages [16 1 (shown in Fig. |5]l considering that only statistical uncertainties are reported and 
that this is not the optimal method for a lifetime measurement. 

To convert the efficiency corrected yields to per event yields requires the true number of NSD 
events, which is obtained by correcting the number of selected events for the event selection 
inefficiency. The event selection includes both the online trigger and offline selection described 
in Section |2j The event selection efficiency is determined in two ways. In the default method, 
it is calculated directly from the Monte Carlo simulation (appropriately weighted by the track 
multiplicity to reproduce the data). In the alternative method, the event selection efficiency ver- 
sus track multiplicity is derived from the Monte Carlo. Then, each measured event is weighted 
by the inverse of the event selection efficiency based on its number of tracks. The number of 
events divided by the number of weighted events gives the event selection efficiency. How- 
ever, since the event selection requires a primary vertex, no events will have fewer than two 
tracks. Therefore, the Monte Carlo is also used to determine the fraction of NSD events which 
have fewer than two tracks and the event selection efficiency is adjusted to include this effect. 
In both methods, the event selection efficiency accounts for unwanted SD events which pass 
the event selection. The numerator in the efficiency ratio contains all selected events, including 
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Figure 5: Kg (left), A (middle), and (right) corrected decay time distributions at v^s = 0.9 and 
7 TeV. The values of the lifetimes, derived from a fit with an exponential function (solid line), 
are shown in the legend along with the world-average value. The error bars and uncertainties 
on the lifetimes refer to the statistical uncertainty only. 

single-diffractive events, while the denominator contains all NSD events. 

5 Systematic uncertainties 

The systematic uncertainties, reported in Table [ij are divided into two categories: normaliza- 
tion uncertainties, which only affect the overall normalization, and point-to-point uncertain- 
ties, which may also affect the shape of the pj and \y \ distributions. 

The list below summarizes the source and evaluation of the point-to-point systematic uncer- 
tainties. 

• Kinematic weighting versus 2D binning: The efficiency corrections using the ID kine- 
matic weighting technique (used for the analysis) and the 2D binning technique 
(used for the V'^ analysis) were compared by measuring the efficiency with both 
methods on the highest statistics channel (Kg at 7 TeV). 

• Non-prompt A: The contribution of non-prompt A decays is varied by a factor of 
two in the simulation. 

• MC tune : The nominal efficiency calculated from the default pythia 6 D6T tune IflSH 
is compared to the efficiency obtained from the pythia 6 PerugiaO (PO) tune |[T7| and 
PYTHIA 8 [18J. 

• Variation of reconstruction cuts : The following cuts are varied for all three modes: 
V*^ vertex separation significance (±2(7"), 3D impact parameter of V^ and (±2c7"), 
3D impact parameter of tracks (±2(7), cut on Kg(A) mass for A(Kg) candidates 
(±1.5(7), and increase of number of hits required on each track from 3 to 5. For the 
E^, additional cuts were varied: the vertex separation significance (±1(7) and E^ 
vertex fit probability (±3%). 

• Detached particle reconstruction: Finding that the corrected lifetime distributions 
are exponential with the correct lifetime is a verification of our understanding of the 
reconstruction efficiency versus decay length. The systematic uncertainty is taken 
as the difference between the fitted lifetimes and the world-average lifetimes I16J . 
While the Kg and A lifetimes are within 1% of the world-average, a 2% systematic 
uncertainty is conservatively assigned. 
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5 Systematic uncertainties 



• Mass fits : As an alternative to using a double-Gaussian signal shape, the V° invari- 
ant mass distributions are fit using a signal shape taken from Monte Carlo. 

• Matching versus fitting: The number of reconstructed events, used in the numerator 
of the efficiency, is calculated in two ways. The truth matching method counts all 
reconstructed candidates which are matched to a generated candidate, based on the 
daughter momentum vectors and the decay vertex. The fitting method fits the MC 
mass distributions to extract a yield. The difference between these two is taken as a 
systematic. 

• Misalignment: The nominal efficiency, obtained using a realistic alignment in the 
MC, is compared to the efficiency from a MC sample with perfect alignment. 

• Beamspot: The location and width of the luminous region of pp collisions (beamspot) 
is varied in the simulation to assess the effect on efficiency. 

• Detector material : The nominal efficiency is compared to the efficiency from a MC 
simulation in which the tracker was modified. The modification consisted of two 
parts. First, the mass of the tracker was increased by 5% which is a conservative 
estimate of the uncertainty. Second, the amounts of the various materials inside 
the tracker were adjusted within estimated uncertainties to obtain the tracker which 
maximized the interaction cross section. Both effects were implemented by changing 
material densities such that the tracker geometry remained the same. The effect is to 
decrease the efficiency as more particles, both primary and secondary, interact. 

• GEANT 4 cross sections : The cross sections used by GEANT 4 for low energy strange 
baryons and all antibaryons are known to be overestimated ||T9) . The size of this 
effect is evaluated by analyzing A-A asymmetries. 

As the trigger efficiency is used to derive the number of NSD events, it only affects the normal- 
ization. The normalization systematic uncertainties, most of which come from trigger efficiency 
uncertainties, are described below. 

• Alternative trigger efficiency calculation: The difference between the default and al- 
ternative trigger efficiency measurements, described in Sec.|4| is taken as the system- 
atic uncertainty on the method. 

• Fraction of SD vs NSD : The change in trigger efficiency when the fraction of single- 
diffractive events in Monte Carlo is varied by ±50% is taken as the systematic un- 
certainty on the fraction of SD events. The pythia 6 MC produces approximately 
20% SD events while the fraction in the triggered data is considerably less ||ll|2l. As 
the UA5 experiment measured 15.5% for this fraction at 900 GeV [.20 J , a variation of 
±50% is conservative. 

• Modelling diffractive events: In addition to the fraction of SD events, the modelling 
of SD and NSD events may not be correct. The trigger efficiency obtained using 
the D6T tune is compared with the trigger efficiency from the PO tune and PYTHIA 
8. In particular, PYTHIA 8 uses a new Pomeron description of diffraction, modelled 
after PHOJET |.2Tll22| , which results in a large increase in the track multiplicity of SD 
events. 

• Track weighting: The track weighting of the Monte Carlo primarily affects the trig- 
ger efficiency. The track weighting requires a measurement of the track multiplicity 
distribution in data and MC. The default track multiplicity distribution is calculated 
from events which pass the trigger, except the primary vertex requirement is not 
applied. Two variations are considered. First, the track multiplicity distribution is 
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measured from events also requiring a primary vertex. As this requires at least two 
tracks per event, the weight for events with fewer than two tracks is taken to be the 
same as the weight for events with two tracks. Second, the track weighting is deter- 
mined with the primary vertex requirement (as in the first case), but without the HF 
trigger. The variation is taken as a systematic uncertainty on the track weighting. 

• Branching fractions: The results are corrected for other decay channels of Kg, A, 
and E^. The branching fraction uncertainty reported by the PDG IIT6l is used as the 
systematic uncertainty. 

The systematic uncertainties at the two centre-of-mass energies are found to be essentially the 
same. The normalization uncertainties and the detached particle reconstruction uncertainty 
are obtained from the average of the results from the two centre-of-mass energies. The other 
point-to-point systematic uncertainties are derived from the higher statistics 7 TeV results. The 
point-to-point systematic uncertainties are measured as functions of pj and \y\ and found to 
be independent of both variables. Therefore, the systematic uncertainties are estimated such 
that they include approximately 68% of the points (representing a la error). The resulting 
systematic uncertainties are summarized in Table [l] 

Table 1: Systematic uncertainties for the Kg, A, and production measurements. 



Source 


K^ (%) 


A (%) 


E- (%) 


Point-to-point systematic uncertainties 








Kinematic weight vs. 2D binning 


1.0 


1.0 


1.0 


Non-prompt A 




3.0 




MC tune 


2.0 


3.0 


4.0 


Reconstruction cuts 


4.0 


5.0 


5.0 


Detached particle reconstruction 


2.0 


2.0 


3.5 


Mass fits 


0.5 


2.0 


2.0 


Matching vs. fitting 


2.0 


3.0 


3.0 


Misalignment 


1.0 


1.0 


1.0 


Beamspot 


1.0 


1.5 


2.0 


Detector material 


2.0 


5.0 


8.0 


GEANT 4 cross sections 


0.0 


5.0 


5.0 


Point-to-point sum 


5.9 


10.7 


12.7 


Normalization systematic uncertainties 








Trigger calculation 


1.8 


1.8 


1.8 


SD fraction 


2.8 


2.8 


2.8 


Diffractive modelling 


1.5 


1.5 


1.5 


Track weighting 


2.0 


2.0 


2.0 


Branching fractions 


0.1 


0.8 


0.8 


Normalization sum 


4.1 


4.2 


4.2 


Overall sum 


7.2 


11.5 


13.4 



For the measurements of dN/dy, dN / dy\y~o, and dN/dpj, the full systematic uncertainty is 
applied. For the A/Kg and S^/A production ratio measurements, the largest point-to-point 
systematic uncertainty of the two particles is used and, among the normalization systematic 
uncertainties, only the branching fraction correction is considered. Note that for the / A 
production ratios, the A branching fraction uncertainty cancels in the ratio. 



10 



6 Results 



6 Results 

The results reported here are normaHzed to NSD inter actioris. The number of NSD raw events 
(given in Sec. |2]) are corrected for the trigger efficiency and the fraction of SD events after the 
selection. The corrected number of NSD events is 9.95 x 10'^ and 37.10 x 10'^ for a/s = 0.9 and 

7 TeV, respectively. 

6.1 Distributions dN/dy and dNIdpi 

The corrected yields of Kg, A, and H^, versus \y\ and pj are plotted in Fig. |6j normalized to 
the number of NSD events. The rapidity distribution is flat at central rapidities with a slight 
decrease at higher rapidities while the px distribution is observed to be rapidly falling. The 
rapidity distributions also show results from three different PYTHIA models: PYTHIA 6.422 with 
the D6T and PO tunes dllZl and PYTHIA 8.135 dl. Fits to the Tsallis function, described 
below, are overlaid on the pj distributions. 

6.2 Analysis of spectra 

The corrected pi spectra are fit to the Tsallis function |i23il , as was done for charged particles IIH 
121. The Tsallis function used is: 

1 dN _ ^ (n-l)(n-2) 
Nnsd dpi nT[nT + m{n - 2)] 

where C is a normalization parameter and T and n are the shape parameters. The results of 
the fits are shown in Table |2] The data points used in the fits include only the statistical uncer- 
tainty. The statistical uncertainties on the fit parameters are obtained from the fit. The system- 
atic uncertainties are obtained by varying the cuts and Monte Carlo conditions (tune, material, 
beamspot, and alignment) in the same way as used to obtain the point-to-point systematic un- 
certainties on the distributions. The systematic uncertainty on the normalization parameter C 
also includes the normalization uncertainty given in Table nl The normalized indicates good 
fits to most of the samples. The T parameter can be associated with the inverse slope parameter 
of an exponential which dominates at low pi, while the n parameter controls the power law 
behaviour at high pj. While both parameters are necessary, they are highly correlated, with 
correlation coefficients around 0.9, making it difficult to elucidate information. Nevertheless, it 
is clear that T increases with particle mass and centre-of-mass energy. This indicates a broader 
low-px shape at higher centre-of-mass energy and for higher mass particles. In contrast, the 
high Pj power-law behaviour seems to show a much steeper fall off for the two baryons than 
for the Kg. While the power-law behaviour of the baryons does not show any dependence on 
the centre-of-mass energy, the fall off of the Kg particles produced at a/s = 0.9 TeV is steeper 
than those produced at a/s = 7 TeV. 

We calculate the average pj directly from the data in the dN / dpj histograms. The Tsallis func- 
tion fit is used to obtain the correct bin centre and to account for events beyond the measured 
Pt range, both of which are small effects. The statistical uncertainty on the average pj is ob- 
tained by finding the standard deviation of px and dividing by the square root of the equivalent 
number of background-free events, where the equivalent number of background-free events is 
given by the square of the inverse of the relative uncertainty on the total number of signal 
events. The systematic uncertainty is composed of two components added in quadrature. The 
first component is the same as used in determining the Tsallis function systematic uncertainties 
(varying the cuts and Monte Carlo conditions). The second component is obtained by using the 



1 + 



Pj + 
nT 



m 



(1) 
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Figure 6: Kg (top), A (middle), and (bottom) production per NSD event versus \y \ (left) and 
Pj (right). The inner vertical error bars (when visible) show the statistical uncertainties, the 
outer the statistical and point-to-point systematic uncertainties summed in quadrature. The 
normalization uncertainty is shown as a band. Three PYTHIA predictions are overlaid on the \y\ 
distributions. The solid curves in the pj distributions are fits to the Tsallis function as described 
in the text. 



12 



6 Results 



Table 2: Results of fitting the Tsallis function to the data. In the C, T, and n columns, the first 
uncertainty is statistical and the second is systematic. The parameter values and ;\;^/NDF are 
obtained from fits to the data with only the statistical uncertainty included. 



Particle 


Vs (TeV) 


C 


T (MeV) 


n 


;^:VNDF 




0.9 


0.776 ± 0.002 ± 0.042 


187±1±4 


7.79 ± 0.07 ± 0.26 


19/21 


A 


0.9 


0.395 ± 0.002 ± 0.041 


216±2±11 


9.3 ±0.2 ±1.1 


32/21 


E" 


0.9 


0.043 ± 0.001 ± 0.006 


250 ± 8 ± 48 


10.1 ±0.9 ±4.7 


19/19 




7 


1.329 ±0.001 ±0.062 


220 ± 1 ± 3 


6.87 ± 0.02 ± 0.09 


50/21 


A 


7 


0.696 ± 0.002 ± 0.058 


292 ± 1 ± 10 


9.3 ±0.1 ±0.5 


128/21 


77 — 


7 


0.080 ± 0.001 ± 0.012 


361 ± 7 ± 72 


11.2 ±0.7 ±4.9 


21/19 



mean pj of the fitted Tsallis function. The average pj from data and PYTHIA 6 with the D6T 
underlying event tune is shown in Table [sj The PYTHIA values are quite close to the a/s = 7 TeV 
data and somewhat lower than the y/s = 0.9 TeV data. Although the average pj results from 
PYTHIA are relatively close to the data, the PYTHIA pj distributions are significantly broader 
than the data distributions. This disagreement can be seen in Fig. |7| which shows the ratio of 
PYTHIA to data for production of Kg, A, and versus transverse momentum. As well as a 
broader distribution, the PYTHIA distributions also show significant variation as a function of 
tune and version. 

Table 3: Average pj in units of MeV/c obtained from the appropriate dN / dpi distribution as 
described in the text. Results from PYTHIA 6 with tune D6T are also given. In each data column, 
the first uncertainty is statistical and the second is systematic. 





v/s = 0.9TeV 


y/s = 7TeV 


Particle 


Data MC (D6T) 


Data MC (D6T) 




654±1±8 580 


790±1±9 757 


A 


837±6± 40 750 


1037±5±63 1071 




971 ±14 ±43 831 


1236 ±11 ±72 1243 



The relative production versus transverse momentum between different species is shown in 
Fig.js] TheN(A)/N(K^) and N(E")/N(A) distributions both increase with pj at low px/as ex- 
pected from the higher average pj for the higher mass particles. At higher pj the N(A) /N(Kg) 
distribution drops off while the N(E^)/N(A) distribution appears to plateau. This is consis- 
tent with the values of the power-law parameter n for these distributions. Interestingly, the 
collision energy has no observable effect on the level or shape of these production ratios. The 
PYTHIA results are superimposed on the same plot. While PYTHIA reproduces the general fea- 
tures, it differs significantly in the details and shows large variations depending on tune and 
version. 

Figure |9] shows a comparison of the CMS px distributions with results from other recent exper- 
iments [3l|2l|25l. To compare with the CMS results, the CDF, ALICE, and STAR distributions 
are multiplied by Snpj, 4, and 87r, respectively. The CDF cross sections are also divided by 
49 mb (the NSD cross section used by CDF [25|) to obtain distributions normalized to NSD 
events, matching the CMS and STAR normalization. The ALICE results are normalized to in- 
elastic events (including single diffractive events). The ALICE and CMS results at 0.9 TeV agree 
for all three particles. The distributions behave as expected, with higher centre-of-mass energy 
corresponding to increased production rates and harder spectra. To remove the effect of nor- 



malization. Fig. 10 shows a comparison of A to Kg and H to A production ratios versus trans 



verse momentum. The CMS results agree with the results from pp collisions at a/s =0.2 TeV 
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Figure 7: Ratio of MC production to data production of Kg (top), A (middle), and (bottom) 
versus px at a/s = 0.9 TeV (open symbols) and \/s = 7 TeV (filled symbols). Results are shown 
for three PYTHIA predictions at each centre-of-mass energy. To reduce clutter, the uncertainty, 
shown as a band, is included for only one of the predictions (D6T) at each energy. This uncer- 
tainty includes the statistical and point-to-point systematic uncertainties added in quadrature 
but does not include the normalization systematic uncertainty. 



14 



6 Results 




Figure 8: N(A)/N(Kg) (left) and N(S-)/N(A) (right) in NSD events versus pj. The inner 
vertical error bars (when visible) show the statistical uncertainties, the outer the statistical and 
all systematic uncertainties summed in quadrature. Results are shown for three pythia pre- 
dictions at each centre-of-mass energy. 



from STAR ["241 and at -\/s =0.9 TeV results from ALICE fS]. These three results show a remark- 
able consistency across a wide variety of collision energies. In contrast, the CDF values for 
N(A)/N(Kg) |26| are significantly higher than the CMS results while the CDF measurements 
of N(E^)/N(A) [25 1 are lower, albeit with less significance. 

Reducing the pj distributions to a single value, the average pj, we compare the CMS results 



with earlier results at lower energies in Fig. 11 13l l24ll26l 432l. The CMS results are in excellent 



agreement with the recent ALICE measurements at 0.9 TeV. The CMS results continue the 
overall trend of increasing average pj with increasing particle mass and increasing centre-of- 
mass energy. 

6.3 Analysis of production rate 

As a measure of the overall production rate in NSD events, ^ lyf^o and the total yield for < 2 



were extracted and tabulated in Table 
the region \y\ < 0.2. The integrated yie 



dy 

The quantity ^lyf^o is the average value of ^ over 
s for lyl < 2 are obtained by integrating the pj spectra. 
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using the Tsallis function fit to account for particles above the measured pj range. 

The central production rates of Kg, A, and are compared to previous results in Fig 
The results show the expected increase in production with centre-of-mass energy with little 
evidence of a difference due to beam particles. As the ALICE results are normalized to all 
inelastic collisions, they are expected to be somewhat lower than the CMS results. 



The production ratios N(Kg)/N(A) and N(E^)/N(A) versus \y\ are shown in Fig. 13 The 
rapidity distributions are very flat and, as observed in the pj distributions of Fig. |8j sEow no 
dependence on centre-of-mass energy. Three pythia predictions at each centre-of-mass energy 



are also shown in Fig. 13 These results confirm what can already be seen in the comparisons 
shown in the left panes of Fig. |6| PYTHIA underestimates the production of strange particles 
and the discrepancy grows with particle mass. 
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Figure 9: Kg (top), A (middle), and (bottom) production per event versus pj. The error bars 
on the CMS results show the combined statistical, point-to-point systematic, and normalization 
systematic uncertainties. The error bars on the CDF p5| , ALICE [3J, and STAR [24J results show 
the combined statistical and systematic uncertainties. The CMS, CDF, and STAR results are 
normalized to NSD events while the ALICE results are normalized to all inelastic events. 



Table 4: ^\y^o arid integrated yields {\y\ < 2.0) per NSD event from data. In each data 
column, the first uncertainty is statistical and the second is systematic. 





V~s = 


0.9 TeV 




7 TeV 


Particle 


dN\ 

dy li/~0 


N 


dN 1 

dy li/~0 


N 




0.205±0.001 ±0.015 


0.784 ±0.002 ±0.056 


0.346 ±0.001 ±0.025 


1.341 ±0.001 ±0.097 


A 


0.108±0.001 ±0.012 


0.404 ±0.004 ±0.046 


0.189±0.001 ±0.022 


0.717±0.005±0.082 


H 


0.011 ±0.001 ±0.001 


0.043 ±0.001 ±0.006 


0.021 ±0.001 ±0.003 


0.080±0.001 ±0.011 
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Figure 10: Ratio of A to Kg production (top) and to A production (bottom) versus pj. The 
CMS, ALICE [3], and STAR |[24l error bars include the statistical and systematic uncertainties. 
The CDF error bars include the statistical uncertainties for N(A)/N(Kg) [26J and the statisti- 
cal and systematic uncertainties for N(E")/N(A) [25J. The CDF N(A)/N(Kg) bin sizes are 
doubled to reduced fluctuations. For experiments in which the binning for A and is differ- 
ent (ALICE and STAR), bins are merged to provide common bin ranges in the N(E^)/N(A) 
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6.3 Analysis of production rate 



17 



> 



o 
> 

CD 
O 



> 

CD 

CD 

w 

I- 
Q. 



0.8 
0.7 
0.6 
0.5 

1.1 
1 

0.9 
0.8 
0.7 
0.6 

1.2 
1 

0.8 
0.6 



Ks 

CMS (pp) 
CDF (pp) 
ALICE (pp) 
UA5 (pp) 
STAR (pp) 



A 

CMS (pp) 
CDF (PP) 
E735 (pp) 
ALICE (pp) 
UA5 (pp) 
STAR (pp) 



CMS (pp) 
E735 (pp) 
ALICE (pp) 
UA5 (pp) 
STAR (pp) 



h 



10" 



CMS 



I 

□ 
I 



I 
I 



1 



I 
I 



s [TeV] 



1 



10 



Figure 11: Average pj for Kg (top), A (middle), and (bottom), as a function of the centre-of- 
mass energy. The CMS measurements are for \y\ < 2. The other results are from UA5 1 
(pp collisions covering \y\ < 2.5, \y\ < 2, and \y\ < 3 for Kg, A, and S^, respectively), E735 1 
(pp collisions using tracks with —0.36 < t] < 1.0), CDF 1261 (pp collisions covering < 1.0), 
STAR [24] (pp collisions covering \y\ < 0.5), and ALICE fSf (pp collisions covering \y\ < 0.75 for 
Kg and A and \y\ < 0.8 for E^). Some points have been slightly offset from the true energy to 
improve visibility. The vertical bars indicate the statistical and systematic uncertainties (when 
available) summed in quadrature. 
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(pp), STAR [24 J (pp), and ALICE [3 J (pp). The CMS, UA5, and STAR results are normalized to 
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offset from the true energy to improve visibility. The vertical bars indicate the statistical uncer- 
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for the CMS, ALICE, and STAR results. 
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Figure 13: The production ratios N(A)/N(K°) (left) and N(H")/N(A) (right) in NSD events 
versus \y\. The inner vertical error bars (when visible) show the statistical uncertainties, the 
outer the statistical and all systematic uncertainties summed in quadrature. Results are shown 
for three PYTHIA predictions at each centre-of-mass energy. 



Table |5] shows a comparison of the production rate of data to pythia 6 with the D6T tune. 
The left column shows a large increase in the strange particle production cross section as the 
centre-of-mass energy increases from 0.9 to 7 TeV. The systematic uncertainties for this ratio 
are reduced as the same uncertainty affects both samples nearly equally. The results for Kg 
and A are consistent with the increase observed in inclusive charged particle production HI 
(ill ~ 1-67) while the results show a slightly greater increase. The increase in particle 
production from 0.9 to 7 TeV is not well modelled by pythia 6. Another feature, seen in the 
right column, is the deficit of strange particles produced by pythia 6. The deficit of Kg particles 
in the MC, 15% (28%) low at 0.9 (7) TeV, is consistent with the results found in the production 
of charged particles However, the deficit is much worse as the mass increases, resulting 

in a 63% reduction in particles in MC compared to data at a/s =7 TeV. While values are 
only presented for PYTHIA 6 with the D6T tune, the same features are also evident for the other 
two PYTHIA comparisons in the rapidity distribution plots in Fig. |6] 



Table 5: Comparison of strangeness production rates between PYTHIA 6 Monte Carlo (D6T) 
and data. In each column, the first uncertainty is statistical and the second is systematic. 



Particle 




\ f ly-o(7TeV) 1 


(D6T) 


v/s = 0.9 


-f |,.o(MCD6T)- 


= 7TeV 


_f |,.o(0.9TeV)_ 
Data MC 


[ f ly-o(Data) J 
TeV 


A 


1.69 ± 0.01 ± 0.06 1.42 
1.75 ± 0.02 ± 0.08 1.48 
1.93 ± 0.10 ± 0.09 1.51 


0.852 ± 0.005 ± 0.061 0.717 ± 0.001 ± 0.052 
0.606 ± 0.007 ± 0.070 0.514 ± 0.003 ± 0.059 
0.477 ± 0.021 ± 0.064 0.373 ± 0.010 ± 0.050 



7 Conclusions 

This article presents a study of the production of Kg, A, and particles in proton-proton col- 
lisions at centre-of-mass energies 0.9 and 7 TeV. By fully exploiting the low-momentum track 
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7 Conclusions 



reconstruction capabilities of CMS, we have measured the transverse-momentum distribution 
of these strange particles down to zero. From this sample of 10 million strange particles, the 
transverse momentum distributions were measured out to 10 GeV/ c for Kg and A and out 
to 6 GeV/c for H~. We fit these distributions with a Tsallis function to obtain information on 
the exponential decay at low pj and the power-law behaviour at high fj. All species show a 
flattening of the exponential decay as the centre-of-mass energy increases. While the baryons 
show little change in the high-pj region, the Kg power-law parameter decreases from 7.8 to 6.9. 
The average px values, calculated directly from the data, are found to increase with particle 
mass and centre-of-mass energy, in agreement with predictions and other experimental results. 
While the PYTHIA pj distributions used in this analysis show significant variation based on 
tune and version, they are all broader than the data distributions. 

We have also measured the production versus rapidity and extracted the value of dN/dy in 
the central rapidity region. The increase in production of strange particles as the centre-of- 
mass energy increases from 0.9 to 7 TeV is approximately consistent with the results for in- 
clusive charged particles. However, as in the inclusive charged particle case, PYTHIA fails to 
match this increase. For Kg production, the discrepancy is similar to what has been found in 
charged particles. However, the deficit between PYTHIA and data is significantly larger for 
the two hyperons at both energies, reaching a factor of three discrepancy for H~ production 
at = 7 TeV. If a quark-gluon plasma or other collective effects were present, we might ex- 
pect an enhancement of double-strange baryons to single-strange baryons and/ or an enhance- 
ment of strange baryons to strange mesons. However, the production ratios N(A)/N(Kg) and 
N(H~ ) /N(A) versus rapidity and transverse momentum show no change with centre-of-mass 
energy. Thus, the deficiency in PYTHIA is likely originating from parameters regulating the 
frequency of strange quarks appearing in colour strings. The variety of measurements pre- 
sented here can be used to tune PYTHIA and other models as well as a baseline to imderstand 
measurements of strangeness production in heavy-ion collisions. 
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